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1.  Introduction 

The  conventional  approach  to  analyzing  antenna  structures  radiating 
in  the  presence  of  a lossy  half -space  involves  repeated  evaluation  of  the 
Sommerfeld  integrals  appearing  in  the  expressions  for  the  vector  potentials 
[1  ].  Since  the  evaluation  of  these  infinite  integrals  is  an  extremely 
time-consuming  process,  much  attention  has  been  focused  in  recent  years  on 
developing  techniques  for  efficiently  evaluating  the  Sommerfeld  integrals 
without  unduly  sacrificing  the  accuracy  [2- 8| . However,  even  the  latest 
reported  techniques  for  evaluating  the  Sommerfeld  integrals  are  40  to  100 
times  slower  than  the  Reflection  Coefficient  Method  (RCM) , which  evaluates 
these  integrals  asymptotically  and  is  valid  for  large  kr  (where  k is  the 
free  space  wave  number,  and  r is  the  distance  between  the  image  and 
observation  points) . 

In  this  paper,  we  present  a new  approach  for  rapid  and  accurate 
numerical  evaluation  of  the  vector  potentials  that  avoids  the  tedious  task 
of  handling  the  Sommerfeld  integrals.  We  begin  with  the  two-dimensional 
Fourier  transforms  of  the  vector  potentials  which  are  conveniently  expressed 
in  simple  closed  forms.  Next,  we  show  that  under  a suitable  approximation 
the  inverse  transform  of  the  vector  potentials  can  be  performed  analytically 
using  a set  of  identities.  The  resulting  space  domain  expressions  are  either 
expressed  in  a closed  form  or  require  evaluating  a finite  integral.  These 
expressions  are  valid  for  a wide  range  of  frequencies,  ground  parameters,  and 
observation  points. 

In  the  following  sections,  the  use  of  the  above  procedure  is  demonstrated 
by  considering  a horizontal  current  element  over  lossy  ground,  and  several 
numerical  examples  are  included  to  illustrate  the  accuracy  and  computational 
efficiency  of  the  method..  We  find  that  the  computational  time  is  only  slightly 
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larger  than  the  reflection  coefficient  method  (RCM)  and  the  accuracy  of  the 
results  is  good  for  a wide  range  of  parameters  of  practical  interest. 

2.  Transformed  Vector  Potentials 

The  fields  radiated  by  a horizontal  current  element  over  a lossy 
ground  (see  Fig.  1)  can  be  expressed  in  terms  of  two  vector  potential 
components  Kx  and  II  [5,8].  We  define  the  following  two-dimensional  Fourier 
transform  pair: 
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{ { n exp [-j (ax  + 3y) ] dx  dy 


(la) 


30 

— j | II  exp  [j  (ax  + 3y)  ] da  d3 
4ir“  -<*> 


(lb) 


where  the  - on  top  represents  the  transformed  quantities.  It  has  been 
shown  [5,8]  that,  using  the  exp(jwt)  time  convention  and  for  observation 
points  above  the  lossy  ground,  the  following  expressions  for  the  vector 
potentials  satisfy  the  Maxwell's  equations  and  the  required  boundary 
conditions : 
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Figure  1.  Geometry  and  the  coordinate  systems  for  the  current 
element  P radiating  over  imperfect  ground,  where 
£ =1  and  £„  * £ - icr/ (<!)£_)  has  been  assumed. 
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xo  = ^Vir*  Idx’ 
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(4b) 
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(4d) 


Equations  (2b)  and  (2c)  are  the  solutions  to  the  problem  of  a current  element 
radiating  in  free  space  and  can  be  interpreted  as  the  direct  and  the 
reflected  contributions  to  the  field  at  the  observation  point  for  a perfectly 
conducting  ground.  Their  inverse  transforms  are  expressed  in  the  following 
well-known  space  domain  forms : 
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(5a) 


nx  " X0  exp(-;j  ki  ri)  / 4ffri 
Kx  " -I0  exp(_J  kl  r2)  / 4ffr2 


(5b) 


where  (r1 , ei , 4*^^)  and  (r2 , 02 , ^>2)  are  the  spherical  coordinate  systems 
erected  at  the  source  and  its  image  point  respectively.  Equations  (2d) 

and  (3)  are  the  correction  terms  to  the  perfect  ground  solution  and  are  not 
directly  amenable  to  inversion  in  closed  forms.  Traditionally,  these  two 

equations  are  transformed  into  space  domain  representations  comprised  of 
infinite  integrals  that  are  oscillatory  in  nature.  These  integral  representa- 
tions take  many  different  forms  and  the  following  are  examples  which 
contain  the  Hankel  functions: 
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The  integrals  appearing  in  equations  (6a)  and  (6b),  better  known  as  the 
Sommerfeld  integrals,  are  quite  time-consuming  to  evaluate  and  represent 
a major  hurdle  to  the  task  of  efficient  analysis  of  antennas  located 
over  a lossy  ground.  Many  of  the  recent  works  in  the  literature  [4-8] 
have  concentrated  on  the  problem  of  reducing  the  computational  time 
required  in  the  evaluation  of  the  Sotnmerfeld  integrals  without  sacrificing 
the  accuracy  beyond  a tolerable  level.  In  the  following  sections,  an 
alternate  approach  is  proposed  for  Fourier  inverting  (2d)  and  (3)  in  a 
more  direct  manner  which  circumvents  the  need  for  computing  the  Sommerfeld 
integrals  (6a)  and  (6b) . 

3.  Approximate  Expressions  for  the  Transformed  Vector  Potentials 

An  examination  of  the  expressions  for  transformed  vector  potentials 
given  in  (2d)  and  (3)  reveals  two  important  and  useful  properties.  First, 
both  of  the  equations  have  an  identical  z-variation  that  corresponds  to  a 
space  domain  solution  emanating  from  the  image  point  of  the  original 
dipole  source.  Second, it  is  apparent  that  the  two  equations  are  well 

behaved  in  the  Fourier  transform  domain  and  decay  exponentially  to  zero 

2 2 2 
outside  the  circle  a +3  ■ k^. 

Starting  from  the  transform  pair  (2c)  and  (5b),  one  can  generate 
a set  of  useful  identities  by  successively  applying  the  operator  3/3z  to 
the  transform  pair.  We  have,  for  instance 

Q * 1 exp [-j  (z+h)] 

Q - 2(j)n+1  g ; n - 0,1,2,... 

-v  n 


(7a) 

(7b) 
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where  g is  che  free-space  Green's  function 


1 


{(*.7.*)  - exp(-J  r2J  / 4irr2  ; r2  - + 7*  + U+h>"  • 


We  also  note  chat  che  successive  partial  derivatives  of  g can  be  obtained 
explicitly;  hence,  (7b)  is  expressible  in  a closed  form.  We  will  now 
attempt  to  approximate  che  transformed  vector  potentials  q.1x  and  ilz> 
given  in  (2d)  and  (3),  such  that  the  inverse  transform  operation  is  performed 
via  (7a-o)  without  an  undue  sacrifice  in  che  accuracy.  To  chis  end,  we 
introduce  che  only  approximation  needed  co  accomplish  chis  goal,  by  meeting 


y2  - « 


Kuo  and  ilai  [7,  Eq.  3]  have  employed  chis  approximation  co  manipulate  che 

Sommerfeld  integrals  (6)  and  have  found  it  co  be  accurate  inside 

o 2 1 

che  i+3  - It’  circle  for  most  practical  parameter  ranges  or  interest. 

In  addition,  out3ide  che  circle  of  visible  range,  che  decaying  exponential 
overcomes  mosc  of  che  errors  introduced  by  che  substitution  of  (9)  into 

(2d)  and  ( 3)  . 


3a.  Approximation  for  ilx 

The  transform  .."x,  given  in  (2d),  can  be  written  in  che  following 


crx  - ro  <vV  Yi  (2+h)1 

(l-<) 


Introducing  the  approximation  in  (9)  and  using  the  identities  in  (7), 


one  can  express  the  space-domain  expression  for  in  a closed  form 

-T  - T 2>k~  3 2 a2 

O'*  0 jkj  (W)  8 ' Iq  k2  a-*)  ,z2* 

Before  closing  this  subsection,  it  may  be  worthwhile  to  point  out 

one  of  the  key  differences  between  the  Kuo  and  Mei  scheme  [7]  and  the  present 

method,  both  of  which  employ  identical  approximations.  In  contrast  to 

the  method  in  [7],  the  present  approach  requires  no  numerical  integration 

and  generates  the  solution  entirely  in  a closed  form. 

U X 


3b.  Approximation  for  Jl7 

The  expression  for  II z , given  in  (3),  can  be  rearranged  into  the 
following  form 

- I0  [!-(<+!)  <Y  YV'y"]  exp [-j  Y;l  (z+h)]  . (12) 

k^<  1 2 

Again  introducing  the  approximation  expressed  in  (9),  one  can  further 
simplify  (12)  into 


IQ  2 exPt-j  Yx  (z+h)]  - 


T c (<+l)  I 

0 .2 
k^  k 


where 


p " exp[-j  vL(z+h) 


c ■ k.//T~ 


(14a) 


(14b) 
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The'  space-domain  expression  for  can  now  be  expressed  in  cerms  of  ?, 
che  inverse  transform  of  ?,  given  in  (14a)  as 


a . t -i.  , . i sisal  ? 

* ° 3 ° ^ < 


Using  che  identities  in  (7),  it  can  be  shown  that  ? satisfies  che  following 
first-order  linear  inhomogeneous  differencial  equation 


21  - jcP  - 21  -OL 

)z  J J 5x  3z 


The  boundary  condition  required  for  the  above  differencial  equation  can  be 
imposed  by  evaluating  che  asympcocic  solution  to  che  vector  oocencial 
expression  in  (15),  derived  in  [5],  at  an  observation  point  which  is 
sufficiently  high  above  ground.  The  asymptotic  solucion  is  given  by 


H -21.  cos  4-  sin  9.  cos  9, 
2a  0 2 2 2 


cos  9?  - *'<-sin,i9. 


<cos  9,  + *<-sin-9., 


• exp(-jk1  r7)/4)rr2 


(17a) 


Using  (17a)  in  (15),  we  arrive  at 


-i  f < 

?U  ’ ■ c(<*l>  i t0  3 


V 3“  (° 


2ir77s<2’> 


where  is  a suitable  height  at  which  the  XC1  approximation  in  (i7a) 
is  valid.  5«'e  point  out  chac  for  che  sake  of  simplicity,  the  x and  y 
dependences  of  the  ?,  " , ar.d  g have  been  suppressed  from  their  arguments. 

The  solution  for  ? can  be  obtained  by  integrating  (16)  and  one  can  arrive 
ac  che  following  form  which  is  convenient  for  numerical  computation  : 
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(z+h) A 

(z'+h) A 

H (present  method) 

2 x 103 

.42 

.50 

1.82  - j.953 

.42 

.75 

1.76  - J1.39 

.42 

1.00 

1.57  - jl.31 

.42 

1.25 

1.64  - jl.22 

.42 

1.50 

1.68  - J1.29 

.42 

1.75 

1.63  - jl.30 

.42 

2.00 

1.62  - jl.25 

At  (z+h) /X  - .42 

3 

Exact  H x JO  ■ 
z 

1.58  - jl.29 

At  (z+h)A  ■ .42 

3 

RCM  H x 10  ■ 

z 

1.64  - j.773 

TABLE  1.  Demonstration  of  the  stability  of  the  present  technique  for 
evaluating  H , as  a function  of  z' . For  this  example, 
e * 10,  a 2 ,01  mhos/m,  frequency  ■ 18  MHz,  r_A  ■ .6, 
and  0 - 45°  . 
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Freq.  . Exact  Sommerfeld  Asymptotic  Evalu-  Present  Method 

(MHZ)  r2/A  Integration  x 10$  at ion  (RCM)  x IQ5  x 1q5 


Freq.  . Exact  Sommerfeld,  Asymptotic  EvaluT  Present  Method 

(MHZ)  2 Integration  x 10  ation  (RCM)  * 1(J  x 10^ 
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P(z)  - (P(z')  - 2j  g (a')]  exp  [jc  (z-z')]  + 2j  ~ g(z) 


-2c  axp(jca) 


— g(z)  exp(-jcz)  dz  . 


(13) 


Several  important  features  of  the  Integral  appearing  in  (13)  will 
now  be  pointed  out.  First,  as  seen  in  Table  1,  the  results  of  integration 
are  quite  stable  as  a function  of  z'.  For  example,  for  z'/X  > 1.0  the 
difference  between  the  exact  H values  and  those  computed  by  using  the 
expression  derived  here  is  less  chan  5%  for  z/X  • 0.42.  Second,,  for  the 
range  of  parameters  investigated  thus  far,  the  integrand  has  been  found  to 
be  quite  smooth.  This  is  due  partly  to  the  fact  chat  the  singularity 
associated  with  g is  located  at  P7  (see  Fig.  1)  and  because  the  range  of 
integration  is  less  chan  one  wavelength.  As  a result  of  these  properties, 

Che  numerical  integration  in  (13)  can  be  carried  out  quite  rapidly. 

Finally,  because  the  parameter  z appears  only  in  the  limit  of  the 
integral,  and  not  in  the  integrand  itself,  the  values  of  vector  potential 
~2  along  an  entire  vertical  line  can  be  rapidly  generated  by  simply  marching 
on  the  incremental  integration  steps.  This  is  in  contrast  to  the  conventional 
integrals  where  the  entire  Integrand  must  be  recomputed  for  each  value  of 
the  observation  point  z . 


4.  Numerical  Results  and  Conclusions 

Tables  2-5  compare  the  accuracy  of  the  present  and  " expressions, 
given  in  (11)  and  (15),  with  those  obtained  via  the  exact  and  asymptotic 
evaluations  of  the  Sommerfeld  integrals.  The  present  technique  is  virtually 
as  accurate  as  the  exact  integration  for  a relatively  high  conducting  ground 
(or  |<!  large)  such  as  sea  water  (see  Tables  2-3)  even  for  extremely  small 
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image  to  observation  point  distance  v^/\.  For  less  conducting  grounds 


(or  | ie  | small),  the  approximation  of  defined  by  (9)  introduces  some 


errors  as  ^A  decreases  (see  Tables  4-5).  However  the  results  are  still 


useful  for  most  practical  antenna  problems  and  remain  superior  to  those 


derived  by  using  the  RCM  method. 


The  procedure  has  been  successfully  tested  for  a wide  variety  of  r^. 


8 , and  |k|  > 5.  On  the  Cyber-175  computer,  the  evaluation  of  a II  and  II 
Z U X z 


pair  for  a given  observation  point  required  * 5 msec  for  the  present  technique. 


while  the  RCM  method  needed  ~ 1 msec,  and  a recent  efficient  Sommerfeld  exact 


integration  technique  [5,8]  typically  required  40-60  msec  of  computing  time. 


Although  not  discussed  here,  it  is  worthwhile  to  mention  that  the 


problem  of  a vertical  dipole  radiating  over  a lossy  ground  can  be  handled 


in  a similar  manner.  The  analysis  of  various  antenna  structures,  comprising 


both  horizontal  and  vertical  wire  sections  over  a lossy  ground,  will  appear 


in  a future  communication. 
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